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Abstract 

It is a well known fact that Dirac phenomenology of nuclear forces predicts the existence of large 
scalar and vector mean fields in matter. To analyse the relativistic self-energy in a model inde- 
pendent way, modern high precision nucleon-nucleon (NN) potentials are mapped on a relativistic 
operator basis using projection techniques. This allows to compare the various potentials at the 
level of covariant amplitudes were a remarkable agreement is found. It allows further to calculate 
the relativistic self-energy in nuclear matter in Hartree-Fock approximation. Independent of the 
choice of the nucleon-nucleon interaction large scalar and vector mean fields of several hundred 
MeV magnitude are generated at tree level. In the framework of chiral EFT these fields are dom- 
inantly generated by contact terms which occur at next-to-leading order in the chiral expansion. 
Consistent with Dirac phenomenology the corresponding low energy constants which generate the 
large fields are closely connected to the spin-orbit interaction in NN scattering. The connection 
to QCD sum rules is discussed as well. 
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I. INTRODUCTION 



A fundamental question in nuclear physics is the role which relativity plays in nuclear 
systems. The ratio of the Fermi momentum over the nucleon mass is about kp/M ~ 0.25 and 
nucleons move with maximally about 1/4 of the velocity of light. This implies only moderate 
corrections from relativistic kinematics in finite nuclei. Non-relativistic approaches such as, 
e.g., Skyrme-Hartree-Fock and relativistic approaches describe finite nuclei equally well. 

However, there exists a fundamental difference between relativistic and non-relativistic 
dynamics: a genuine feature of relativistic nuclear dynamics is the appearance of large scalar 
and vector mean fields, each of a magnitude of several hundred MeV. The scalar field S s 
is attractive and the vector field is repulsive. In relativistic mean field (RMF) theory, 
both, the sign and the size of the large scalar and vector fields are enforced by the nuclear 
saturation mechanism At nuclear saturation density po — 0.16 fm~ 3 the empirical fields 
deduced from RMF fits to finite nuclei are of the order of E s ~ —350 MeV and S — +300 
MeV (In mean field theory only the time-like component of contributes in static 
systems with time-reversal symmetry). 

A problem is, however, that these scalar/vector fields are no direct observables as, e.g., 
the nuclear binding energy or the nucleon potential. The single particle potential in which 
the nucleons move originates from the cancellation of the two contributions U cent ~ S + S s 
and is of the order of -50 MeV. Therefore one has no direct experimental access to the 
interpolating scalar/vector fields. There exist, however, several features in nuclear structure 
which can be explained naturally within Dirac phenomenology while models based on non- 
relativistic dynamics have difficulties or, at least, one has to introduce additional model 
parameters. The most well known feature is the large spin- orbit splitting in finite nuclei. In 
a relativistic framework the strong spin-orbit force appears naturally from the coupling to 
the lower components of the Dirac equation where the scalar-vector mean fields add up in 
the spin-orbit potential Us.o. oc (S — S s ) — 750 MeV. Due to this fact RMF theory is able 
to reproduce the strong spin-orbit splitting in spherical nuclei quantitatively without the 



introduction of additional parameters [2|. A second symmetry, observed more than thirty 
years ago in single-particle levels of spherical nuclei is the so called pseudo-spin symmetry [3] . 
While all attempts to understand this symmetry within non-relativistic approaches failed, it 
can naturally be understood within RMF theory as has been shown by Ginocchio J] a few 
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years ago. This symmetry, again a consequence of the coupling to the lower components, is 
exact in the limit So = — S s and is broken in nature by the amount (Eo + £ s ) / (So — S s ) which 
is less than 10%. A third example are the moments of inertia in rotating nuclei. Relativistic 
dynamics implies that in the rotating system a Coriolis term occurs due to the spatial vector 
currents, however, with all couplings already fixed through the time-like components Q]. 

An alternative approach for nuclear matter are ab initio many-body calculations. Based 
on high precision nucleon-nucleon (NN) interactions one treats short-range and many-body 
correlations explicitly. A typical example for a successful many-body approach is Brueckner 
theory . In the relativistic Brueckner approach the nucleon inside the medium is dressed 
by the self-energy X. The in-medium T-matrix is obtained from the relativistic Bethe- 
Salpeter (BS) equation and plays the role of an effective two-body interaction which contains 
all short-range and many-body correlations of the ladder approximation. Solving the BS 
equation the Pauli principle is respected and intermediate scattering states are projected out 



of the Fermi sea. The summation of the T-matrix over the occupiec 
sea yields finally the self-energy in Hartree-Fock approximation ^ 



states inside the Fermi 
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contrast to relativistic Dirac-Brueckner-Hartree-Fock (DBHF) calculations which came up in 
the late 80ties non- relativistic BHF theory has already half a century's history 1(3] . Despite 
strong efforts invested in the development of improved solution techniques for the Bethe- 
Goldstone (BG) equation, the non-relativistic counterpart of the BS equation, it turned out 
that, although such calculations were able to describe the nuclear saturation mechanism 
qualitatively, they failed quantitatively. Systematic studies for a large number of nucleon- 
nucleon interactions showed that saturation points were always allocated on a so-called 
Coester-line in the E/A — p plane which does not meet the empirical region of saturation. 
In particular modern one-boson-exchange (OBE) potentials lead to strong over-binding and 
to too large saturation densities where relativistic calculations work much better [lH llq. 

However, in relativistic approaches the nuclear interaction is always described in some 
sort of a meson exchange picture. The mesons represent effective bosonic degrees of freedom 
which are either directly adjusted to the properties of nuclear matter and finite nuclei, as in 
the case of RMF theory, or to vacuum NN scattering. Hence it is a fundamental question 
to decide whether the large scalar and vector fields enforced by Dirac phenomenology of 
nuclear systems are an artefact of the meson exchange picture or whether they reflect a 
deeper characteristics of nature. 
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A connection to Quantum-Chromo-Dynamics (QCD) as the fundamental theory of strong 
interactions is established by QCD sum rules (f| Q|. The change of the chiral condensates 
(qq), {q^q) in matter leads to attractive scalar and repulsive vector self-energies which are 
astonishingly close to the empirical values derived from RMF fits to the nuclear chart. 

It is remarkable that relativistic many-body calculations yield again scalar and vector 
fields which are of the same sign and magnitude as obtained from RMF theory or, alter- 
natively, from QCD sum rules. Such a coincidence could not have been expected a priori. 
Moreover, DBHF calculations jl4l agree even on a quantitative level surprisingly well with 
the QCD based approach of Ref. [2| where chiral fluctuations from the long and intermediate 
range pion-nucleon dynamics were considered on top of the chiral condensates. 

These facts suggest that preconditions for the existence of large fields in matter or, alter- 
natively, the density dependence of the QCD condensates, must already be inherent in the 
vacuum AW interaction. The connection of the nucleon-nucleon force to QCD is given by 
the fact that the interaction is described by the exchange of the low lying mesonic degrees 
of freedom. The long-range part of the interaction is mediated by the one-pion-exchange 
(OPE) while the scalar isoscalar intermediate range attraction is mainly due to correlated 
two-pion-exchange. The short-range part, i.e., the hard core, is dominated by light vec- 
tor meson exchange, i.e., the vector isoscalar uo meson and the vector isovector p. Modern 
one-boson-exchange potentials (OBEP) as e.g. the Bonn potentials [17[ are based on the 
exchange of these mesons and provide high precision fits to nucleon-nucleon scattering data. 
Meson-nucleon coupling constants and form factors are empirically fixed from the data. 
Thus OBEPs are the result of relativistic phenomenology at the level of the elementary AW 
interaction. There exist, however, also high precision non-relativistic empirical potentials 

n n 

such as the Argonne potential |18[ or the Nijmegen potentials |2Cj. 

A more systematic and direct connection to QCD is provided by chiral effective field 
theory (EFT). Up to now the two-nucleon system has been considered at next-to-next-to- 

nnn 

next-to-leading order (N 3 LO) in chiral perturbation theory |22j, I23L l24| . In such approaches 
the AW potential consists of one-, two- and three-pion exchanges and contact interactions 
which account for the short-range contributions. The advantage of such approaches is the 
systematic expansion of the AW interaction in terms of chiral power counting. The ex- 
pansion is performed in powers of (Q/A x ) u where Q is the generic low momentum scale 
given by the nucleon three-momentum, or the four-momenta of virtual pions or a pion mass. 
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A x ~ Anf^ ~ 1 GeV is the chiral symmetry breaking scale which coincides roughly with 
the Borel mass A# of the sum rules. In such an expansion the low-energy constants (LECs) 
related to pion-nucleon vertices can be fixed from pion-nucleon scattering data j^J. 

A better understanding of the common features and the differences of the various ap- 
proaches is essential in order to arrive at a more model independent understanding of the AW 
interaction, in particular since all the well established interactions fit AW scattering data 
with approximately the same precision. A direct comparison of relativistic phenomenology 
based on the meson exchange picture with chiral EFT and non-relativistic phenomenology 
is, however, difficult since the latter two approaches lack of a clear Lorentz structure. At low 
momentum scales the different potentials can be mapped on each other using renormaliza- 
tion group methods |25j]. This led recently to the construction of a "model independent" low 
momentum potential V\ ow k by integrating out the dynamics for momenta above a cut-off 
scale of about A ~ 2 fm _1 25]. It has been argued that beyond this scale the short-range 
part of the interaction, mediated by vector meson exchange or point-like counter terms, 
becomes dominant and leads to the deviations of the various approaches. 

Although a breakthrough in some sense, the renormalization group approach does not 
help to clarify the relativistic structure of the potentials which is essential e.g. in order to 
generate (or not to generate) large scalar/vector mean fields in nuclear matter. 

The present work tries to answer this question. We apply projection techniques to map the 
various potentials on Dirac phenomenology. The philosophy behind this approach is based 
on the fact, that any AW interaction, independent whether relativistic or non-relativistic, 
contains a certain spin-isospin operator structure. By projection techniques this operator 
basis is mapped on the operator basis of Dirac phenomenology which is given by the Clifford 
algebra in Dirac space. This allows to identify the different Lorentz components of the 
interaction. Starting from the angular-momentum representation of a given AW potential, 
one transforms to plane-wave helicity states and finally to Lorentz invariant amplitudes in 
Dirac space (3, I23]. Such a transformation is well defined in the positive energy sector 
for on-shell amplitudes and allows to compare the AW potentials on the basis of Lorentz 
invariant amplitudes. A remarkable agreement between relativistic and non-relativistic OBE 
potentials, non-relativistic phenomenological potentials and EFT potentials, respectively, 
has been found in 4f|. This agreement is also reflected in the structure of the relativistic 
self-energy when we further calculate the mean field in infinite nuclear matter in Hartree- 
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FIG. 1: (Color online) Diagonal matrix elements V(q, q) in the 1 Sq partial wave for different high 
precision NN potential models. 



Fock approximation at tree-level. The scalar and vector self-energy components are found 
to be large, i.e., of the order of several hundred MeV 4||. The present work extends the 
investigations of Ref. The formalism is outlined in detail and we discuss the density 

dependence of the fields as well as the implications for the nuclear equation-of-state. The 
connection between chiral EFT and QCD sum rules is investigated. The present formalism 
allows a quantitative extraction of the scalar /vector fields which are generated by pion 
dynamics and contact terms at different orders in the chiral expansion. 

The paper is organised as follows: in sec. II we discuss the operator structure of the 
various potentials. The transformation onto the covariant basis is outlined in sec. Ill where 
also the results of this analysis namely the Lorentz invariant amplitudes are shown. Sec. 
IV contains the determination of the relativistic mean fields in nuclear matter. In sec. V 
the structure of the relativistic self-energy fields from chiral EFT is discussed, as well as the 
connection to QCD sum rule predictions. Finally the self-consistent Hartree-Fock results for 
the equation of state calculated with three different potentials (Bonn A, Nijm93 and Nijm I) 
are discussed in sec. VI. 
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II. OPERATOR STRUCTURE OF THE NN POTENTIALS 



A. OBE potentials 



As typical examples for modern high precision OBEPs we consider the Bonn A 37] and 
the high-precision, charge-dependent Bonn potential (CD-Bonn) |28j]. The Bonn potentials 
are based on the exchange of the six non-strange bosons (tt, i], p, u, S, a) with masses below 
1 GeV. These are the two scalar mesons a (isoscalar) and S (isovector), the two pseudo- 
scalar mesons ir (isovector) and rj (isoscalar), and the two vector mesons u (isoscalar) and 
p (isovector). The potentials are derived in the no sea approximation which neglects the 
coupling to anti-particles. 

The Born scattering matrix is given by the sum over the corresponding scalar, pseudo- 
scalar and vector mesons and has the following structure 

vtf, q) = E i) 4 2) D a (q> - q) /#> . (1) 

a=s,ps,v 

In the two-nucleon centre-of-mass frame (cm.) the four- momenta of the incoming nucleons 
are q9~' 2 ' = (E(q), ±q) and correspondingly, the four-momenta of the outgoing nucleons are 
g^ 1 / 2 ) = (.E'(q'), ±q'). The initial and final relative cm. momenta are g M = \{qy — qy) 
and q' = \{q'f' — q'^), respectively. For on-shell scattering |q| = |q'| with E(q) = -E(q') = 
\/ M' 2 + q 2 the energy-transfer is zero, i.e., q'^ — q^ = (0, q' — q). The matrices (0) factorise for 
each meson a into the form factors T a at each meson-nucleon vertex, the meson pro pag ator 
D a and the meson-nucleon vertices K a themselves. In the standard Bonn potentials [37j the 
phenomenological form factors have the form 

/ A 2 - m 2 \ na 

where m a is the corresponding meson mass and A a is a cut-off in order to avoid divergences 
at short distances. The meson propagators read 

D (q'-q)-i I D *>U q) ± ^ = *™ Z *1M (3 ) 

for scalar and pseudo-scalar mesons s,ps and vector mesons v. The Dirac structure of the 
potential is contained in the meson-nucleon vertices 

Ks ~ (2vr) 21 ' Kps ~ (2vr) 2 2M %1 ' (2vr) 2 ^ 7 + 2M m ) ■ {A} 



For the pseudo-scalar mesons n and r\ a pseudo-vector coupling is used in order to fulfil 
soft pion theorems. The vertices of the isovector bosons n, 5, p obtain additional r 2 ■ T\ 
isospin matrices which are suppressed in Eqs. (JU). The uj meson has no tensor coupling, i.e., 

= 0. 

The potential, i.e., the OBE Feynman amplitudes are obtained by sandwiching V between 
the incoming and outgoing Dirac spinors 

V(q',q)= £ ^(q',q) D a (c^ - q) U 2 (-^K^ 2 (-q) U^)^^^ • (5) 

ct=s,ps,v 

The relativistic operator structure is thus completely determined by the matrix elements of 
the vertices k q . In helicity representation the Dirac spinor basis is given by 



E + M 1 . 

MA(q) = V^iW~ juisL |xa ' (6) 

V E+M 

where \\ denotes a two-component Pauli spinor with A = ±|. The normalisation of the 
Dirac spinor is chosen such that U\U\ = 1. 

A consequence of the Feynman amplitudes (0) is their general non-local structure which 
distinguishes the field theoretical relativistic OBE approach from local non-relativistic po- 
tentials. This is even true for the relativistic OPE compared to the local, non-relativistic 
OPE (see e.g. the discussion in j^). However, for on-shell scattering the relativistic am- 
plitudes acquire a local structure in the sense that they are functions of q 2 and q' — q. 
In particular for forward and backward scattering, i.e., 9 = 0,tt, the amplitudes are "local" 
functions of q 2 and q. The non-local structure of the relativistic amplitudes becomes evident 
when going off-shell, e.g. in the intermediate states in the Bethe-Salpeter equation 

n 

The standard Bonn (A,B,C) potentials |3j| contain 13 free parameters for coupling con- 
stants and cut-off masses and two additional parameters if one considers the masses of the 
scalar mesons as effective parameters. The matrix elements are calculated with the OBNNS 
code of R. Machleidt [3^ when Bonn A is used and the corresponding CDBONN package 
of R. Machleidt when CD-Bonn is used. 

In contrast to the standard Bonn potentials [37| the OPE part of the CD-Bonn poten- 
tial accounts for charge symmetry breaking in nn, pp and np scattering due to the 
different pion masses m n o and m^±. The CD-Bonn potential can be referred to as a phe- 
nomenological NN potential, since by fine-tuning of the partial wave fits \ 2 P er datum is 
minimised to 1.02, adding up to a total of 43 free parameters. 
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1. Non-relativistic reduction 



The OBE potentials as e.g. the Bonn potentials can be reduced to a non-relativistic 
representation by expanding the full field-theoretical OBE Feynman amplitudes into a set 
of spin and isospin operators 

V = Y^i + ylrvr 2 \O i . (7) 

i 

The operators Oj obtained in this low energy expansion, assuming identical particle scatter- 
ing and charge independence, are defined as 

01 = l, 

02 = f i • cr 2 , 

3 = (ovk)(<r 2 .k) (8) 

4 = |(<Tl + cr 2 ) ■ n, 

5 = (cr r n)(<x 2 -n), 

where k = q' — q, n = qxq' = Pxk and P = |(q + q') is the average momentum. 
The potential forms Vi are then functions of k, P, n and the energy. In order to perform a 
non-relativistic reduction, usually the energy E is expanded in k 2 and P 2 

/k 2 \* k 2 P 2 

E(q) = — + P 2 + M 2 ~M+ — + — . (9) 

and terms to leading order in k 2 /M 2 and P 2 /M 2 are taken into account. The meson propa- 
gators D a (k 2 ) given in Eq. (jSJ) are approximated by their static form (—1)/ (k 2 + m 2 ). The 
equivalent to Eq. (JZJ) in configuration space is given by 

01 = l, 

2 = c r cr 2 , 

03 = #12 = 3((Ti-r)(o- 2 -r) - <Tl-(T2, (10) 

4 = L • S, 

5 = Q12 = |[(tr r L)(o- 2 -L) + (<T 2 -L)(<rrL)]. 

These operators are the well known central, spin-spin, tensor, spin-orbit and quadratic spin- 
orbit operators, respectively. The total angular momentum is denoted by L = r x P and 
the total spin S = \{cr\ + cr 2 ). 



B. Non-relativistic potentials 



1. Meson-theoretical potentials 
We consider the modern Nijmegen soft-core potential Nijm93 [19j] as the first example o 
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a non-relativistic meson-theoretical potential. It is an updated version of the Nijm78 
potential, where the low energy NN interaction is based on Regge-pole theory leading to the 
well known OBE forces. The contributions considered in this model are the pseudo-scalar 
mesons 7r, 77, rj', the vector mesons p, <p, uo and the scalar mesons 5, S*, e and the Pomeron P 
and the J = tensor contributions, leading all in all to a number of 13 free parameters. Since 
it is constructed from approximate OBE amplitudes it is based on the operator structure 
given in Eq. (jHJ) plus an additional operator Oq = \(cri — <x 2 ) • L accounting for charge 
independence breaking which is new compared to the older version Nijm78. Exponential 
form factors are used. This potential gives a x 2 P er datum of 1.87, which is comparable to 
similar OBE potentials like the standard Bonn potentials. 



2. Phenomenological potentials 

Another class of non-relativistic NN potentials are the so called high quality poten- 
tials where x 2 l^data ~ 10. Here we study the Nijmegen potentials Nijm I, Nijm II and 
Reid93 Q|. The Nijm I and Nijm II potentials are both based on the Nijm78 potential. In 
the Nijm I potential some nonlocal terms in the central force are kept whereas in the Nijm II 
potential all nonlocal terms are removed. Although based on the meson-theoretical Nijm78 
potential these potentials are often referred to as purely phenomenological models, since the 
parameters are adjusted separately in each partial wave leading to a total of 41 parameters. 
At very short distances, both potentials are regularised by an exponential form factor. 

The Nijmegen soft-core Reid93 [3] potential is a phenomenological potential and is 
therefore based on a completely different approach. In the meson-theoretical Nijmegen 
potential Nijm93 the potential forms Vi are the same for all partial waves, whereas in the 
Reid93 potentials every partial wave is parametrised separately by a convenient choice of 
combinations of central, tensor and spin-orbit functions (local Yukawas of multiples of the 
pion mass) and the related operators, i.e., the operators 0\ to O4 from Eq. (fTU|) . It is 
regularised by a dipole form factor and has 50 phenomenological parameters giving all in all 
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a x 2 /Ndata = 1-03. All the Nijmegen potentials contain the proper charge dependent OPE 
accounting for charge symmetry breaking in nn, pp and np scattering due to different pion 
masses m w o , m n ± . 

The same holds for the Argonne potential v\$ ^J, also an example for a widely used 
modern high precision phenomenological NN potential. It is given by the sum of an electro- 
magnetic (EM) part, the proper OPE, and a phenomenological intermediate- and short-range 
part unrestricted by a meson-theoretical picture: 

V = V EM + V" + V R . (11) 

The EM interaction is the same as that used in the Nijmegen partial- wave analysis. Short- 
range terms and finite-size effects are taken into account as well |l8j |. 

The strong interaction part V n + V R can be written in a form like given in Eq. (jZJ) in 
configuration space, where the Argonne v±s potential is not constructed by approximating 
the field-theoretical OBE amplitudes (except for the OPE), but by assuming a very general 
two-body potential constrained by certain symmetries. The potential forms V{ parametrising 
the intermediate and short-range part are mostly local Woods-Saxon functions. 

The local two-body operators are the same charge independent ones used in the Argonne 
v 14 potential 

O t = l,a v a 2 ,S l2 ,L-S,L 2 ,L 2 {a v a 2 ),(L-S) 2 . (12) 

Due to isovector exchange these operators have to be multiplied by the isospin matrices 
7~i ■ r 2 which than adds up to 14 operators. Additionally, four operators accounting for 
charge independence breaking are introduced 

0i=15,18 = ^12) (Cl-C2)^12) S 12 T U , [j zl + T z2 ) , (13) 

where T\ 2 = 3t z1 t z2 — Ti-t 2 , is the isotensor operator, defined analogously to the spin tensor 
S12 operator. 

Thus the operator structure is more general than that imposed by a non-relativistic, local 
OBE picture, in particular for the intermediate and short distance part. In total, Argonne 
t>i8 contains 40 adjustable parameters and gives a y 2 per datum of 1.09 for 4301 pp and np 
data in the range 0-350 MeV . The code used to calculate the potential matrix elements 
of the Argonne fig model in momentum space was provided by H. Muether and T. Frick. 
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C. Low-Energy potentials 



1. EFT potentials 

Following the concept originally proposed by Weinberg 32] there has been substantial 
progress in recent time in order to derive quantitative NN potentials from chiral effective 
field theory. As already mentioned, the chiral expansion is performed in powers of (Q/A x ) v 
where v — corresponds to leading order (LO), v — 2 to next-to- leading order (NLO), v = 3 
to next-to-next-to-leading (N 2 LO) and finally v = 4 to next-to-next-to-next-to-leading order 
(N 3 LO). It turned out that for a quantitative description of NN scattering data one has 



to go up to N 3 LO 22J, |23|, |24j in the chiral expansion for the two-nucleon problem. N 2 LO 
contributions were still found to be very large compared to NLO. This implies that 2ir (and 
37r) contributions have to be included up to order four. The effective chiral Lagrangian can 
be written as 

r - r(2) i i r( 2 ) i r( 3 ) . hzn 

where the superscript refers to the number of derivatives or pion mass insertions (chiral 
dimension) and the ellipsis stands for terms of chiral order four or higher. The corresponding 
chiral NN potential is then defined by 

{sum of irreducible 
\ + contacts. (15) 
7T + 2tz contributions I 

The 27r exchange contributions to the iViV interaction at order four have been derived by 
Kaiser j^j]. Recently, quantitative iViV potentials including contact terms at order four were 
derived by Entem and Machleidt, the so-called Idaho potential 22, and by Epelbaum, 
Glockle and Meissner 2^ |. 

For the present comparison we apply the Idaho potential [23]. The operator structure of 
the momentum-space NN amplitude has the general form given in Eq. (jZJ) with the operators 
Oi from Eq. (jSJ). The potential forms Vi (i = C, S, T, LS, aL) can be expressed as functions 
of |(q' — q)| and |k|. 

The Idaho potential is regularised by an exponential cut-off 

y(q', q ) V (q', q) e^'/^ e^ 2 " (16) 
where A = 0.5 GeV in all partial waves. This does not affect the chiral order of the potential, 
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but introduces contributions bey ond that order. The total number of free model parameters 
in the N 3 LO potential is 29 [23J. 

For the evaluation of the matrix elements we applied the N 3 LO program package provided 
by D. Entem and R. Machleidt. 

2. Renormalization Group approach to NN interaction 

Recently, another approach has been proposed to arrive at a better model independent 
understanding of the NN interaction 25j] . In this approach a low-momentum potential 
Viow k is derived from a given realistic NN potential by separating the low-momentum part, 
i.e., by integrating out the high-momentum modes, and using renormalization group (RG) 
methods to evolve the NN potential models from the full Hilbert space to the low momentum 
subspace. At a cutoff of A = 2.1 fm _1 all the various NN potential models were found to 
collapse to a model-independent effective interaction Vj ow k- 

Since elastic NN scattering data constrains the NN interaction only up to a momentum 
scale of about 400 MeV, which corresponds to the pion threshold, modern high precision 
potentials differ essentially in the treatment of the short-range part, as depicted in Fig. 
The philosophy behind the RG approach is to replace the unresolved short distance structure 
by something simpler, e.g. contact terms, without distorting low-energy observables. 

III. TRANSFORMATION TO A CO VARIANT OPERATOR BASIS 
A. Covariant operators in Dirac space 

Any two-body amplitude can be represented covariantly by Dirac operators and Lorentz 
invariant amplitudes. A detailed discussion of the general structure of relativistic two-body 



amplitudes can be found in Refs. 
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34j . However, a relativistic treatment invokes automati- 
cally the excitation of anti-nucleons. Nucleon-nucleon scattering, in both, the non-relativistic 
approaches discussed above but also in the framework of the standard OBE potentials is 
restricted to the positive energy sector and neglects the coupling to anti-nucleons. As a 
consequence one has to work in a subspace of the full Dirac space which leads to on-shell 
ambiguities which require some care. 



13 



The inclusion of negative energy excitations with 4 states for each spinor yields altogether 
4 4 = 256 types of two-body matrix elements with respect to their spinor structure. Symme- 
try arguments reduce these to 44 for on-shell particles 34|. If one takes only the subspace 
of positive energy solutions into account this leads to 2 4 = 16 two-body matrix elements. 
Considering in addition only on-shell matrix elements the number of independent matrix 
elements can be further reduced by symmetry arguments down to 5. Thus, all on-shell two- 
body matrix elements can be expanded into five Lorentz invariants. These five invariants are 
not unique since the Dirac matrices involve always also negative energy states. Therefore 
a decomposition of the one-body NN potential into a Lorentz scalar and a Lorentz vector 
contribution depends to some part on the choice of these five Lorentz invariants. 

A natural choice of a set of five linearly independent covariant operators to represent 
a 4 x 4 Dirac matrix are the scalar, vector, tensor, axial-vector and pseudo-scalar Fermi 
covariants 



1(8)1, V = 7 m ® 7m , 



T = a^ u <x> a 



/IV ; 



7s® 7s, A = 757^^757^. (17) 



Since one works with physical, i.e., antisymmetrized matrix elements, one has to realize 
that the Fierz transformation T [26] couples direct and exchange covariants which mixes 
the different Lorentz structures 
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(18) 



The covariants on the left hand side of Eq. (ITH1) are the interchanged Fermi covariants 
defined in Ref. |26J as 



S = SS, V = SV, T = ST, A = SA, P = SP, 



(19) 



where the operator S exchanges the Dirac indices of particles 1 and 2, i.e., Sm(1) (7 m(2) t = 
u(l) T u(2) a . Therefore the direct covariants T m with m = {S,V, T,P,A} can be expressed 
in terms of the exchange covariants r m with m = {S,V,T,P,A}. 

In contrast to the NN potentials where the pion-nucleon coupling is given by a pseudo- 
vector vertex, the set (jl7H9J) contains the pseudo-scalar covariant P. This suggests to replace 
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P in Eqs. (|1 7119(1 by the corresponding pseudo-vector covariant 

pv = ^i 75 ® ti 75 . (20) 

2M /5 2M r5 V ; 

This leads to an on-shell equivalence since the matrix elements of the pseudo-vector and 
the pseudo-scalar matrix operators are identical in the case of on-shell scattering between 
positive energy states: 

^(q , )^75«(q)=«(q')75n(q) . (21) 

On the other hand the PV vertex suppresses a coupling to antiparticles since the overlap 
matrix elements vanish for on-shell scattering 

^(q')^75^(q) = . (22) 

In order to identify the PV contributions clearly in the antisymmetrized amplitudes - note 
that due to the Fierz transformation (II 811 all operators are coupled - one can switch to a set 

n 

of covariants originally proposed by Tjon and Wallace |34j |. Based on the following operator 
identities 

-(T + T) = S + S + P + P (23) 
V + V = S + S- P- P (24) 
one finds that the following set of covariants 

r m = {S,S,(A-A),PV,PV} (25) 



provides a set of Dirac operators for the positive energy sector |34j which completely sepa- 
rates the direct and exchange pv contributions from the remaining operator structure. This 
has the advantage that the OPE exchange which is dominant at low energies is decoupled 
from the remaining amplitudes and gives only a contribution to the PV operator. In the 
following we will refer to the set of covariants in Eq. (|23jl as the pseudo — vector representa- 
tion and that of Eq. (|17jl as the pseudo — scalar representation. Note that on-shell matrix 
elements of PV, PV in (125(1 are equivalent to those where the pseudo-vector covariants are 
replaced by P, P. 

The on-shell equivalence does not affect physical observables which are built on complete 
matrix elements as e.g. the single particle potential U 

C/(k) s . p . ex £ (u(k)ii(q)|V(k,q)|u(k)u(q) - u(q)u(k)) (26) 
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but it leads to uncertainties in operators which are, like the self-energy X, based on traces 
over only one particle. As discussed in |35| . a pseudo- vector nN coupling leads to the 
pseudo — vector representation ()25|1 as the most natural choice of the relativistic operator 
basis. 



37| |. The last step depends on the 



B. Projection onto the covariant operators 

In this section the technique is described necessary to project the Born amplitudes from 
an angular-momentum basis onto the covariant basis, given by Eqs. (j!7|) or (|25|). The 
procedure is standard and runs over the following steps 

\LSJ) — > partial wave helicity states — > plane wave helicity states — > covariant basis. 

The first two transformation can be found in Refs. |3f 

choice of the covariant operator basis, see e.g. [3, I27J]. Here we sketch the essential steps 
briefly. 

Independent of the various models, the amplitudes are determined normally in the 
|LS'JM)-representation and can be denoted as V^^(q',q). In case of on-shell scattering 
(|q| = |q'|), due to time- reversal invariance and spin and parity conservation, only five of 
sixteen possible matrix elements are linearly independent for a fixed total angular momen- 
tum J (spin singlet and triplet states). By inversion of Eq. (3.32) in these five partial 
wave amplitudes are transformed from the |LS'JM)-representation into the partial wave 
helicity representation | JMA1A2) and are then decoupled via inversion of Eq. (3.28) from 
Ref. 36]. Since we deal with two-nucleon states which are two-fermion states, we have to 
evaluate the fully antisymmetrized matrix elements by restoring the total isospin I = 0, 1 
via the standard selection rule 

= _i. (27) 

The five plane wave helicity matrix elements are then obtained by a summation over the 
total angular momentum J 

(A^q'l^l AxA 2 q) = £ (^^) ^ A ,(0)<A;Ayv^V,q)|AiA 2 > . (28) 

Here 9 denotes the scattering angle between q' and q while A = Ai — A2 and A' = X[ — A' 2 
denote the in- and outgoing helicity states. The reduced rotation matrices d J xx ,{6) are those 
defined by Rose [38|. 
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These plane wave helicity matrix elements can now be projected onto a set of five covariant 
amplitudes in Dirac space. A set of five linearly independent covariants is sufficient for such 
a representation since on-shell we deal with five matrix elements independent of the chosen 
representation. Using the covariants of Eq. (|17|) (the 'pseudo-scalar choice') the on-shell 
potential matrix elements for definite isospin / can be represented covariantly as 27 1 

^(|q|,60 = Fl(\q\,0)S + F*(\q\,0)V + F*(\q\,0)T 

+F*(|q|,0)P + Fl(|q|,0)A . (29) 

The Lorentz invariant amplitudes i^(|q|, 9) with m = {S, V, T, P, A} from Eq. (|2"9"|) depend 
only on the relative cm. momentum |q| and the scattering angle 9 and are related to the 
plane wave helicity states defined in Eq. ()28|) by 

(A' 1 A^q / |y I |A 1 A 2 q) = ^(A' 1 A' 2 q / |r m |A 1 A 2 q)Fi(|q|^) . (30) 

m 

The indices (1) and (2) refer to particle one and two. Eq. (J5UJ) is a matrix relation between 
the five independent plane wave helicity amplitudes V} (where i = {X[, A' 2 , Ai, A 2 } = 1, 5 
denotes five of sixteen possible amplitudes) and the five unknown covariant amplitudes 
F^(|q|,#). For fixed values of the variables (|q| = |q'|,#) this equation can be written in a 
more compact form 

Vi = TxE c *»fi- ( 31 ) 

The covariant amplitudes are obtained by matrix inversion of Eq. (|31|) which corresponds 



to Eq. (3.23) of Ref. |27|. 

Eq. (fBTj) has to be inverted for two scattering angles, i.e., for 9 = for the direct and 
9 = 7r for the exchange part of the interaction. These two scattering angles are required 
for the Hartree-Fock potential. Details of the inversion of Eq. (}3*Tj) . as well as the treat- 
ment of kinematical singularities of the matrix Cj m occurring at 9 = and 9 = it are given 
in appendix C of Ref. [27. where Eq. (|3*Tj) is explicitly given for 9 = and 9 = ix (Eqs. 
(CIO, 11)). Following Ref. 27( we calculate the real part of the five Lorentz invariant ampli- 
tudes F^ o,1 (|q|,0 = 0,7r) for the direct and exchange case in both, the isospin singlet and 
triplet channels. When derived from physical partial wave amplitudes which are already 
antisymmetrized according to the selection rule (|27|). the exchange amplitudes F m (|q|,7r) 
contain redundant information. 
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Since we are restricted to the subspace of positive energy states, the choice of a set of 
five linearly independent covariants suffers from on-shell a mbig uities, as discussed above. 



Thus the set of covariants 



scattering matrix reads 



is a more appropriate choice |l4|. In this representation the 



V\\q\,9) = ^(|q|^)S-^(|q|^)S + ^(|q|^)(A-A) 
+^v(|q|^)PV-^(|q|,e)PV . 



(32) 



The new amplitudes g l m are related to the Lorentz invariant amplitudes from Eq. 
by the linear transformation 

( 9l \ 
i 
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(33) 



As mentioned before, the representation of the potential given in Eq. (|3*2*|) has the advan- 
tage that the OPE contribution to the amplitudes is completely decoupled from the rest of 
the interaction. The OPE contributes only in the pseudo-vector exchange amplitude Q—^ 1 
and vanishes in all other amplitudes g<$ FE = 5 , g )PE = 9a PE = 9pv E = 0- Thus one avoids that 
the low momentum behaviour of these four amplitudes is to large extent dominated by OPE 
exchange contributions which are present in all five amplitudes from Eq. (j2Uj) due to the 
Fierz transformation. In order to compare the various potentials at the level of covariant 
amplitudes the pseudo-vector representation is therefore the most efficient and transparent 
one. 



C. Covariant amplitudes 

In order to demonstrate the dependence of the relativistic amplitudes on the choice of the 
operator basis we consider in Fig. |2]first the single OPE. The figure shows the corresponding 
amplitudes F m of the pseudo-scalar representation (fT7j) and the g m amplitudes of pseudo- 
vector representation (|23|) . both for the OPE part of the Bonn A potential. Since we 
are dealing with antisymmetrized amplitudes it is sufficient to consider the direct Lorentz 
invariants F m (|q|, 9 = 0) and <? m (|q|, 9 = 0) at scattering angle 9 = 0. As the starting point 
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FIG. 2: Covariant amplitudes for the single OPE for the different choices of the relativistic operator 
basis, i.e., the pseudo-scalar representation F m (\q\,9 = 0) (left) and the pseudo-vector representa- 
tion g m (\q\,6 = 0) (right). 



the OPE is given in the \LSJ) basis and antisymmetrization is ensured by the selection rule 
(|27j). The figure shows the isospin averaged amplitudes denned as 

F m (|q|,0) := \ [i£=°(|q|,0) + 3i£ =1 (|q|, 0)] (34) 

and correspondingly for g m . It is evident that in the pseudo-scalar representation all ampli- 
tudes F m have large non-vanishing contributions from OPE due to the mixing of direct and 
exchange contributions described by the Fierz transformation (|18|). Moreover, as discussed 
above the on-shell equivalence for the pseudo-scalar covariant P and the pseudo-vector co- 
variant PV in leads to identical Lorentz invariant amplitudes Fps = -Fpv = Fp Q|- 
The pseudo- vector representation ([25)1. on the other hand, has the advantage that it decou- 
ples the OPE contribution from the remaining amplitudes, i.e., the OPE gives a non-zero 
contribution only in the amplitude while the others are zero. For the single pion ex- 
change g~ is now easy to interpret: it is just the pion propagator times the pion-nucleon 
form factor (|2*|). 

When the various NN potentials are compared, this is done most efficiently in the pseudo- 
vector representation. All potentials contain an OPE of similar strength which dominates 
at small momenta. The pseudo-vector representation decouples the OPE contribution from 
the remaining amplitudes g m ^ g~ and allows thus a more transparent investigation of the 
short and intermediate range parts of the potentials which are actually the interesting ones. 
Fig-Elshows the isospin-averaged amplitudes <7^(|p|, 9 = 0) for Bonn A, CD-Bonn, Argonne 
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FIG. 3: (Color online) Isospin-averaged Lorentz invariant amplitudes g^(|q|, = 0) for the different 
NN potentials after projection on the Dirac operator structure. The pseudo-vector representation 
of the relativistic operator basis is used. As a reference the amplitudes from solely OPE and from 
a + to exchange, both with Bonn A parameters, are shown. 
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vis, Nijm93, Nijmegen I and II, Reid93, the effective low momentum interaction V\ ow k and 
the chiral Idaho potential. The amplitudes are obtained going through the transformation 
scheme discussed above. Partial waves are taken into account up to J = 90 (Bonn A, 
CD-Bonn, Idaho), J = 9 (Argonne t>i 8 , Nijmegen I/II, Nijm93, Reid93) and J = 6 (V[ ow k)- 

The amplitudes determined from the complete NN potentials are no more as easy to 
interpret as for a single meson exchange where they represent essentially the propagators 
times the form factors. This is also true for the full OBE since the contributions from the 
various mesons are coupled through their exchange parts. Since these amplitudes are not 
very transparent quantities, Fig. El includes as a reference in addition the contributions from 
only OPE and from only a and u exchange, both taken from Bonn A. 

Several features can now be seen from Fig. EJ First of all the four amplitudes g$, g$, 
and g~ are very close for the OBEPs Bonn A, CD-Bonn and Nijm93 and the phenomeno- 
logical non-relativistic Argonne Vis and Nijmegen I/II potentials. Only at very small |q| 
Argonne vis shows a deviating structure. The direct pseudo-vector amplitude gpy falls 
somewhat out of systematics. This amplitude is, however, of minor importance since it does 
not contribute to the Hartree-Fock self-energy (|41H43j) and to the single particle potential. 

The dominance of the OPE at low |q| is reflected in the pseudo- vector exchange amplitude 
g~ which is at small |q| almost two orders of magnitude larger than the other amplitudes. 
In the OBEPs the high momentum part of the interaction, on the other hand, is dominated 
by heavy meson exchange and the corresponding amplitudes gs, 9& approach the a + uj 
exchange result. Deviations from the a + uj amplitudes, e.g. due to exchange of isovec- 
tor mesons p and 5 in the OBEPs are moderate at large |q|. These deviations are more 
pronounced at small |q|. 

The remarkable agreement between the OBE amplitudes and those derived from the 
non-relativistic Argonne Vi$ potential demonstrates two things: first of all, it means that 
for on-shell scattering the Argonne v is can be mapped on the relativistic operator structure 
where the local phenomenological functions V{, Eq. (JZj), play the same role as the meson 
propagators plus corresponding form factors in the meson exchange picture. Secondly, the 
effective treatment of the short- distance physics in Argonne v is is very similar to that in 
the OBE potentials Bonn A, CD-Bonn and Nijm93. This fact can be estimated from Fig. ^ 
where the 1 Sq partial wave amplitudes are close as well. On the other hand the softer 
character of the Reid93 and also the Nijmegen I and II potentials is reflected clearly in the 
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stronger deviation from the a + 00 amplitudes at large |q|. 

Finally we are turning to the effective low momentum potentials V| ow k and the chiral 
Idaho N 3 LO potential. V\ ow k is only shown up to the intrinsic cut-off of 400 MeV. In this 
momentum range the amplitudes fall practically on top of those from the Idaho N 3 LO po- 
tential. At low |q| the amplitudes derived from Idaho N 3 LO and Vj ow k behave qualitatively 
and quantitatively like the previous ones, i.e., they are very close to Bonn A, CD-Bonn and 
Argonne v\$. We conclude that also the effective low momentum potentials can be mapped 
on a relativistic operator structure. For the Idaho N 3 LO potential which is also based on 
the operator structure given in Eq. (jHJ), the functions V\ and V( in combination with the 
corresponding operators, derived from fourth order 2n exchange plus contact terms, lead to 
a structure which is similar to that imposed by the OBE picture. However, clear deviations 
appear in the cut-off region between 400 and 500 MeV. The short-range interactions are 
strongly suppressed by the exponential cut-off form factors and as a consequence the Idaho 
approaches rapidly the OPE result for momenta above 400 MeV. 



IV. SELF-ENERGY IN NUCLEAR MATTER 



With the covariant amplitudes at hand, one is able to determine the relativistic mean 
field in nuclear matter with its scalar and vector components. To do so, we calculate the 
relativistic self-energy £ in Hartree-Fock approximation at tree level. We are thereby not 
aiming for a realistic description of nuclear matter saturation properties which would require 
a self-consistent scheme. Moreover, short-range correlations require to base such calculations 
on the in-medium T-matrix rather than the bare potential V. This leads to the relativistic 
Dirac-Brueckner-Hartree-Fock scheme which has been proven to describe nuclear saturation 
with saving accorac, QQHQ. The seif-consistent iteration of the 

self-energy in combination with the Dyson equation for the in-medium nucleon propagator 
and the Bethe-Salpeter equation for the in-medium T-matrix leads to self-energy components 
which are qualitatively of similar magnitude than the tree level results, as will be seen later 
on. 

The self-energy is determined by the summation of the interaction of a nucleon with 
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four-momentum k with all nucleons inside the Fermi sea in Hartree-Fock approximation 



S a/ g(fc, k F ) 



G?„(q) [v(| P |,o) aCT;/3r -nipu; 



(27T) 



(35) 



Since we work with fully antisymmetrized matrix elements which contain already the direct 
(Hartree) and exchange (Fock) contributions, it is sufficient to evaluate the Hartree integral 
for the self-energy 



E a p(k, k F ) 



lhj g?m [V A (\p\,0) aa ., 



(2tt 



(36) 



G (q) is the Dirac propagator describing the on-shell propagation of a nucleon with mo- 



mentum q inside the Fermi sea in the nuclear matter rest frame 



G u (q) = [i + M}2m5(q 2 - M 2 )Q(q )Q(k F - |q|) 



(37) 



The G functions account for the fact that only positive energies are considered. Here, k, 
taken along the z-axis, is the single particle momentum of the incoming nucleon in the nuclear 
matter rest frame. The relative momentum in the two-nucleon cm. frame where the matrix 
elements V are evaluated, is given by |p| = \J s/A — M 2 , where s = (E(k) + E(q)) 2 — (k + q) 2 
is the total energy of the two nucleons. 

Using the pseudo- vector representation for the on-shell matrix elements V, Eq. fl32(l . the 
self-energy operator reads 



E a p(k,k F ) 



rf 3 q Q(k F - |q|) 
(2vr) 3 4E7(q) 



(ft a 



4M 2 



+ml a/3 



4gs - 9§ + 4 9 'a 



(k<* - q") 



^2 



AM 2 



-98 + 2 9a 



(Jfe" - q 11 ) 2 
AM 2 



-9 



PV 



(38) 



Translational and rotational invariance, hermiticity, parity conservation, and time reversal 
invariance determine the Dirac structure of the self-energy j^. In the nuclear matter rest 
frame the self-energy can be written as 



E(k, k F ) = S s (fc, k F ) - 7o E (/c, k F ) + 7 • k£ v (/c, k F ) 



(39) 



Note that the sign convention for the vector field E = S s — 7^^^ with = (So, kS v ) in 

, Q • It differs from that used standardly 



Eq. (J39)) is that used standardly in DBHF 
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in QHD (£ = S s + 7 M £ M ) and also that of Eqs. (f4T?j) and (foTj). The self-energy components 
are Lorentz scalar functions depending on the Lorentz invariants k 2 , k ■ j and j 2 , where 
denotes the four- vector baryon current. In nuclear matter at rest the time-like component 
is just the baryon density and spatial components of the current vanish, i.e., j M = (p#, 0). 
Hence, the Lorentz invariants can be expressed in terms of k , |k| and kp, where kp denotes 
the Fermi momentum. The components of the self-energy are computed by taking the 
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39] 



respective traces in the Dirac space 
S s = htr [£] , S = 



In doing so, the Lorentz components of the self-energy operator (JH8j) are given by 



— tr [ 7o E] , S v = j^pfr [7 • k£] ■ 



(40) 



S s (/c, k} 



and 



£ v (/c, k} 



d 3 q 

(2^)3 
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(41) 



(42) 



Qz 



AM 2 



PV 



(43) 



(2tt) 3 v " '^'Ikl^q) 

In Fig. 0] the tree level scalar and vector self-energy components in nuclear matter are 
shown obtained with the various iViV potentials at nuclear saturation density with Fermi 
momentum k-p = 1.35 fm^ 1 which corresponds to a density of p = 0.166 fm -3 . As a 
remarkable result, all potentials yield scalar and vector mean fields S s and Eo of comparable 
strength: a large and attractive scalar field E s ~ —(450 -j- 400) MeV and a repulsive vector 
field of — Eo — +(350 4- 400) MeV. These values are comparable to those derived from RMF 
phenomenologically and also from QCD sum rules. Also the explicit momentum dependence 
of the self-energy is similar for the various potentials. The Idaho mean fields follow the other 
approaches at low k but show a stronger decrease above k ~ 2 fm -1 which reflects again 
the influence of the cut-off parameter. Fig. El shows the spatial component of the vector 

re various potentials agree quite well. As known 



self-energy k£ v , Eq. (j43j) . Also here t 
from self-consistent DBHF calculations 



12 



141 ] . the spatial vector self-energy is a moderate 



correction to the large scalar and time-like vector components E s and E . This is found to be 
also the case at tree level where k£ v is about one order of magnitude smaller than the other 
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FIG. 4: (Color online) Tree level scalar and vector self-energy components in nuclear matter at 
kp = 1.35 fm _1 obtained with different NN interaction models. 
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FIG. 6: (Color online) Density dependence of the tree level scalar and vector self-energy components 
in nuclear matter obtained with the various potentials. 



two components. The spatial self-energy originates exclusively from exchange contributions, 
i.e., the Fock term, and vanishes e.g. in the mean field approximation of RMF theory. 

Fig. El displays the density dependence of the fields, evaluated at momentum k = kp. At 
moderate densities the different potentials yield scalar and vector fields which are rather 
close in magnitude. At higher densities the results start to split up which reflects again the 
different treatment of short distance physics in the various interactions. Only the two low 
momentum interactions Idaho N 3 LO and V\ ow k lie practically on top of each other. In this 
context we want to stress again that these results are obtained in lowest order in density. 
Hence, the results are only 'realistic' in the low density limit but not at higher densities 
since short-range correlations are missing. 

In order to estimate the influence of short-range correlations and self-consistency, in Fig. [7| 
the tree level result from Fig. 0] for Bonn A to a corresponding full DBHF calculation are 
compared at kp = 1.35 fm _1 . For DBHF the approach of ^| is used (subtracted T-matrix 
in pv representation). The DBHF calculation yields reasonable saturation properties with a 
binding energy of -Ebind = —15.72 MeV and a saturation density of p = 0.181 fm~ 3 It is 
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FIG. 7: (Color online) Tree level scalar and vector self-energy components in nuclear matter 
at kp = 1.35 fm _1 are compared to corresponding values from a full self-consistent relativistic 
Brueckner (DBHF) calculation. In both cases the Bonn A potential is used. 

no doubt that higher order correlations are essential for saturation of nuclear matter. The 
correlations lead to a general reduction of the vector self-energy by a shift of about 70 MeV. 
Self-consistency and correlations also weakens the momentum dependence, in particular for 
E s . However, except of the 70 MeV shift of S , the absolute magnitude of the self-energies 
is not strongly modified in the realistic calculation. This means that one can expect that 
the large attractive scalar and repulsive vector mean fields will also persist for the other 
interactions when short-range correlations are accounted for in a full relativistic many-body 
calculation. 

Fig. |H] shows finally the single particle potential in nuclear matter at kp = 1.35 fm~ , 
determined from the relativistic self-energy components. The single particle potential is 
defined as the expectation value of the self-energy 



U s . p .(k, k F ) 



< u(k)\j°J:\u(k) > 
< u{k)\u{k) > 



E(k) 



M 



< u(k)\L\u(k) > 



(44) 
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and reads 

f/ s D (k, fc F ) = — L s ^— = — . — S H — ; • (45) 

pA E E Vk 2 + M 2 Vk 2 + M 2 V ; 

Eq. (}4"5j) represents the single particle potential at tree level, i.e., the expectation value of E 
with the bare spinor basis. The next step towards a self-consistent treatment would be to 
use an in-medium spinor basis which includes the scalar and vector self-energy components 
via effective masses and effective four-momenta 

M*(k,k F ) = M + E a (k,k F ), k; = k^ + J:^k,k F ) . (46) 

This would, however, involve higher order corrections in the baryon density and is not 
intended in the present investigations which are restricted to leading oder. 

The single particle potential reflects the well known fact that phase-shift equivalent two- 
body potentials which describe NN scattering data with about the same accuracy [30], 
can be rather different [30j. This can already be seen from Fig. ^ where the 1 5'o matrix 
elements of the various potentials are shown. The differences are mainly due to a different 
treatment of the short-range part of the nuclear interaction, i.e., the hard core which is 
not well constraint by scattering data. Thus the various potentials lead to about the same 
T-matrices when iterated in the Lippmann-Schwinger or Bethe-Salpeter equation. However, 
at tree-level the hard core contributes fully to U s , p , which explains the shift of the various 
results in Fig. |HJ Integrating out the high momentum components, e.g. by renormalization 
group methods, one arrives at equivalent low-momentum potentials V\ ow k 3- Since V\ ow k 
contains no significant contributions from the hard core it gives already at tree level a realistic 
single-particle potential. The situation is similar for the chiral EFT N 3 LO Idaho potential. 
As can be seen from Fig. [T] Idaho is rather close to Vj ow k, not only in the 1 5'o partial wave, 
and correspondingly both lead to comparable potentials. However, the slight shift of about 
10 MeV between Vj ow k and Idaho reflects again the subtle cancellation effects between the 
large scalar/vector fields, since at the scale of the fields, Fig. |3 both lie practically on top 
of each other. 

In the present context the single particle potential serves as an important check of the 
whole procedure. In Fig. |H1 the single particle potential U s _ p , is shown, calculated from 
Eq. (j45j) . i.e., after projecting the NN potentials from the partial wave basis onto the 
covariant operator basis, determining then the relativistic self-energy components and fi- 
nally C/ S .p.. Fig. |H1 includes also the results from a 'non-relativistic' calculation of Z7 B . p , where 
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FIG. 8: (Color online) Single particle potential in nuclear matter at k-p = 1.35 fm _1 , determined 
from the tree level Born amplitudes of the various potentials. The single particle potential deter- 
mined from the relativistic self-energy components after projection onto the covariant operator basis 
is compared to a non-relativistic calculation (stars) where partial wave amplitudes are summed up 
directly. 



the partial wave amplitudes are directly summed up. To do so we used a non-relativistic 
Brueckner-Hartree-Fock program j^] and determined the single particle potential in Born 
approximation. The non-relativistic results are represented by stars in Fig. |S] and shown up 
to a momentum of 400 MeV. This avoids distortions from non-relativistic kinematics which 
occur at higher momenta. At moderate momenta the non-relativistic and the relativistic 
calculations show an excellent agreement which demonstrates the accuracy of the applied 
projection techniques. One has thereby to keep in mind that U s , p , originates in the relativis- 
tic approach from the cancellation of the two scalar and vector fields which are both of the 
order of about 400 MeV. 
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V. THE STRUCTURE OF THE SELF-ENERGY FROM CHIRAL EFT 



With the projection formalism at hand one is now able to investigate the connection 
between the appearance of the matter fields and chiral dynamics in more detail. It allows in 
particular a straightforward and transparent discussion of the contributions which arise at 
different orders in the chiral expansion of the NN interaction, see Eqs. (fT4"Jl and (|15p. Such 
an investigation allows also to build the bridge to the reduction of the in-medium quark 
condensates which is usually interpreted as a signature for a partial restoration of chiral 
symmetry. 



A. Role of contact terms 



We are now in the situation to calculate the relativistic scalar and vector self-energies from 
a chiral EFT nucleon-nucleon potential order by order. For this purpose we apply again the 
chiral Idaho potential [44]. This allows to separate the contributions from different orders 
in the chiral expansion of the NN interaction and provides a connection to the low energy 
constants (LECs) which appear at the different orders. 

Fig. El shows the tree level results for the scalar and vector self-energy components in 
nuclear matter at kp = 1.35 fm _1 obtained in leading order (LO) up to next-to-next-to- 
next-to-leading order (N 3 LO). 

To leading order the chiral NN interaction does not generate significant mean fields. 
The scalar self-energy S s is of the order of about -70 MeV and the vector self-energy is 
practically zero. At LO only the static OPE and contact terms without derivatives appear 
which involve the operators 0\ and O2 from the operator basis (JHJ). Hence at LO no pieces 
from vector exchange occur which would involve all operators Oi, i = 1..5. The small scalar 
field means, on the other hand, that the nucleon mass M*, Eq. I[46j) . does not change 
significantly in matter to leading order in chiral EFT. The dominant contributions arise at 
next-to-leading order (NLO). NLO involves leading two-pion-exchange (2PE) and contact 
terms with two derivatives. The NLO contact terms contain the full operator structure 0{. 
At this level both, scalar and vector self-energy components of about =f400 MeV magnitude 
are generated. Also the signs, i.e., the attractive scalar and the repulsive vector mean field, 
are fixed at NLO. The higher orders, N 2 LO and N 3 LO provide corrections which tend to 
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FIG. 9: (Color online) Tree level scalar and vector self-energy 



1.35 fm 1 obtained with the chiral EFT NN interaction 



components in nuclear matter at 



23j. The fields obtained in leading 



order (LO) up to next-to-next-to-next-to-leading order (N 3 LO) are shown. 



reduce the NLO result, are, however, moderate. N 2 LO contains subleading 2PE and no 
contact terms at all, while N 3 LO contains sub-subleading 2PE, leading three-pion-exchange, 
corrections to OPE and 2PE and contact terms with four derivatives 24 1 . 



In order to investigate the role of pion dynamics and that of contact terms in more detail, 
Table |l] contains the contributions which arise from pion dynamics £ w , i.e., OPE, 2PE, 
3PE and corrections, and those from the contact terms E^ cont ) separately. The contributions 
to the self-energy at a particular order is given by the sum + £( cont ), the full self- 
energy at a certain order v is obtained by adding the contributions from the lower orders 
Yjy) — Y?\=q51 ■ From Table |l] it becomes evident that the dominant contributions to 
the scalar and vector self-energy are generated by the contact terms which arise at next- 
to-leading order. At N 2 LO no contact terms occur in the chiral expansion. The N 3 LO 
contacts provide sizeable corrections to both, scalar and vector self-energy components and 
are of opposite sign than the NLO contributions. The contribution from pion dynamics to 
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S s 




v (cont) 


-S 


_Zj o 


v (cont) 
_2j 


LO 


-64.76 


17.14 


-81.9 


4.49 


19.02 


-14.53 


NLO 


-344.22 


4.4 


-348.62 


376.47 


5.16 


371.31 


N 2 LO 


2.06 


2.06 





-41.92 


-41.92 





N 3 LO 


56.82 


-89.34 


146.16 


-43.27 


79.06 


-122.33 


sum 


-350.1 


-65.74 


-284.36 


295.77 


61.32 


234.45 



TABLE I: Contributions from pion dynamics and contact terms to the scalar and vector self-energy 
components (in MeV) which appear at different orders in the chiral expansion. The evaluation is 
performed at nuclear saturation density hp = 1.35 fm - . 

the self-energy components are found to be generally moderate. The largest contributions 
appear at N 3 LO and are of opposite sign than those from corresponding contact terms. 

Hence the reduction of the nucleon mass M* = M+E s is driven by short- distance physics, 
dominantly by contact terms which occur at NLO. These are four-nucleon contacts with two 
derivatives. At this order the short-range spin-orbit interaction (proportional to O4 in (JBJ) 

zC 5 (o- 1 + cr 2 )-(qxq / ) (47) 

is generated. The appearance of large scalar/vector fields at NLO is therefore in perfect 
agreement with Dirac phenomenology where the large spin-orbit force is intimately connected 
to the appearance of the scalar /vector fields which are generated by short-range isoscalar 
scalar (a) and vector meson (a?) exchange 14, I^J. In EFT the strength of the short- 



range spin-orbit interaction is determined by the C5 parameter which is given by a linear 



combination of the 3 P-wave low energy constants (LECs) 23J, |24 ] 



1 



16tt 



[2C*3po + 3C3P1 — 5C3P2] 



(48) 



Hence the short-range spin-orbit interaction is dictated by P-wave NN scattering. As shown 
by Kaiser the large values of the C5 parameter is in good agreement with corresponding 
values extracted from high precision OBE type potentials (Bonn, CD-Bonn, Nijm93, Ni- 
jmegen I, II) and from Argonne Vis which are all in the range of 3C5/8 ~ 80 4- 90 MeV fm 5 . 
In |£2| these values were also compared to purely phenomenological Skyrme type density 
functionals designed for nuclear structure calculations The values of the correspond- 

ing spin-orbit strength parameter Wq in Skyrme models are also very close in magnitude, 
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FIG. 10: (Color online) Influence of the C5 low energy constant. The figure compares the self- 
energies at NLO to those where all contacts except of C5 are switched off and those results where 
the strength of the C5 parameter is varied. 



i.e. 3Wo/4 ~ 75-^97 MeV fm 5 . The contribution from chiral OPE to the spin-orbit terms in 
the density functional were found to be almost negligible (less than 1%). The lowest order 
irreducible 2PE which occurs at NLO in the chiral expansion provides moderate corrections 
to the iso-scalar spin-orbit strength function whereas the iso-vector strength is more strongly 
affected (2PE contributions lead to a ~ 30% reduction) [52]. Thus the analysis of Kaiser is 
fully consistent with the small fields and £q of ~ =f5 MeV generated by pion dynamics 
at NLO, as observed within the framework of the present analysis. 

Fig. [1JI1 analyses the dependence of the fields on the value of the C5 low energy constant 
in more detail. As already mentioned, at LO two contact terms (C% and C2) appear and at 
NLO, respectively, 5 contacts (C3 to C7). The figure contains the full NLO result, including 
contributions from LO and NLO pion dynamics and contacts and compares this to the case 
where all contacts which appear up to NLO were switched off except of the C5 contribution. 
It contains in addition results with again all contributions, however, scaling the value of C5 
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down to 50%, 10% and 0.1%. It becomes evident that the large scalar and vector mean 
fields are a direct consequence of the large value of C5. Chiral EFT is therefore not only 
in qualitative but quantitative agreement with the picture known from meson-exchange. In 
both cases the fields are related to short distance physics and their strength is dictated by 
P-wave NN scattering data where the spin-orbit forces occur. 



B. Connection to QCD sum rules 

In finite density QCD sum rules scalar and vector fields arise naturally from the structure 
of the quark propagator which is proportional to the corresponding condensates. As shown 
by Cohen et al. [6( the quark correlation function can be expressed to leading order in terms 
of the scalar condensate (p\qq\p) already present in vacuum, and the vector condensate 
(p\q q\p) which is introduced by the breaking of Lorentz invariance due to the presence 
of the medium. The identification of the correlation function with the in-medium nucleon 
propagator of a dressed quasi-particle leads to scalar and vector self-energies S s and E 
which are of the same order in the condensates 

^ 87r 2 r , . . . . ., 8ir 2 m ... 

s s = —ptVspMp) - (??>] = -tt^ttztps ( 49 ) 



- s o = - y^r(p\qioq\p) = — -j^-p ■ (so) 



64vr 2 . ._ 1 . 32tt 2 

— {pmoq\p) = — tt 

B B 

These expression are of leading order in density, ps in (J4l^) is the scalar nucleon density, 
f n = 93 MeV the weak pion decay constant and m u ^ are the current quark masses of 
about 5 4- 10 MeV. The pion-nucleon sigma term = (N\m u uu + rriddd\N) is determined 
by the u and d-quark content of the nucleon and represents the contribution from explicit 
chiral symmetry breaking to the nucleon mass through the small, but non- vanishing current 
quark masses. It has an empirical value of about <jn — 50 MeV. The Borel mass scale 
Ab — 4:7i f n ~ 1 GeV is the generic low energy scale of QCD which separates the non- 
perturbative from the perturbative regime. It coincides with the chiral symmetry breaking 

n 

brmula js| for the nucleon mass M ~ —jrr(qq) one 



scale A x of ChPT. Applying Ioffe's 
finally obtains the fields in the form 



01 



S s (p) = -^ri^Ps , (51) 



A(m u + m d )M 

" Eo(p) = <n " ■ (52) 
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However, the dependence of the nucleon mass in matter on the quark condensate is not as 
straightforward as expression (JoTj) suggests. Concerning the in-medium condensate one has 



careful 



origin 



y 



to distinguish between contributions from the pion cloud and those of non-pionic 
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As pointed out by Birse |41| a naive direct dependence of the nucleon mass on the 
quark condensate through Eq. ()51|) leads to contradictions with chiral power counting. The 
contributions from low momentum virtual pions which enter the in-medium condensate 
should not contribute by the same amount to the change of the nucleon properties in matter. 
They can therefore not as easily be associated with a partial restoration of chiral symmetry 
as the mean field field approximation, Eqs. (J5T1 152j) . would suggest. This problem has also 
been investigated by Chanfray et al. [42] in the framework of the linear sigma model. In 
their studies the authors were able to reconcile the phenomenology of Quantum Hadron 
Dynamics with chiral theory, in that case the linear sigma model. Their conclusion was 
that, in contrast to the scalar condensate (p\qq\p) which is driven by the sigma field, i.e., 
the chiral partner of the pion, the lowering of the nucleon mass M* is driven by a chiral 
invariant scalar field which corresponds to fluctuation along the chiral circle. With other 
words, the condensate is to large extent reduced by the pion cloud surrounding the nucleons 
while the nucleon mass is not. 

To set up the context for the following discussion, we shortly sketch the argumentation of 
Birse j^]J: From Eq. (JBTj) follows that the effective nucleon mass M* = M + S s (p) is directly 
proportional to the nucleon sigma term 

M * = M ( 1 -^) • (83) 

The chiral expansion of the sigma term leads to |43j 

A 2 9 ( 97TNN\ 2 3 

a N = Am v - — — — m v + ... (54) 



16vr V 2M 

In the chiral limit the pion-nucleon coupling is connected to the axial vector coupling by 
the Goldberger-Treiman relation g n NN = 9aM/ f n . The coefficient A involves counter terms 
related to short- distance physics whereas the non-analytic 0{m\) term arises purely from 
long-distance physics of the pion cloud. Inserting ([5*4"J) into (JHHJl implies a dependence of the 
effective nucleon mass M* on the pion mass which is of order 0(m n ). 

At the mean field level, i.e., in T — p approximation, the scalar self-energy (}4*Tj) is on the 
other hand given by the scalar forward scattering amplitude T s (q = 0) ( T s (q = 0) in (J55J) 
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corresponds to the direct amplitudes Fs and gs in and (JH^j). respectively.) 

S s (fc F ) = T s (q = 0)p . (55) 

A comparison of Eq. (|B3j) with Eqs. (|B^j) and (|5^|) would imply that the scalar part of the 
forward scattering amplitude contains a constant and a term of order m^. Such a dependence 
contradicts, however, chiral power counting. In chiral EFT the leading term in the pion mass 
in the NN interaction originates from the low energy expansion of the OPE and is of order 
0(m^) [2^. I23I \l\ . Hence the NN interaction cannot contain a term directly proportional 
to a N /fl 

For the comparison of the sum rule predictions we turn to the density dependence of the 
self-energy. Fig. ITT1 shows the density dependence of the fields from the various orders. As in 
Fig- El the scalar S s , time-like vector E and spatial vector S v self-energies are determined 
at momentum k — kp. The density dependence is shown up to k-p = 1.8 fm _1 which 
corresponds to about 2.5 times nuclear saturation density. As can be seen from Fig. ITT1 the 
relative contributions from the various orders remain the same over the entire density range 
considered. For comparison the figure contains also the corresponding fields as predicted 
by leading order QCD sum rules, i.e., Eqs. (joTjl and (|5^|) . For the evaluation of Eqs. (J5Tj) 
the empirical value of = 50 MeV has been chosen for the nucleon sigma term, f n = 93 
MeV and (m u + m^) = 12 MeV. For the evaluation of the scalar field in (|51jl we have set 
the scalar density equal to the vector density, i.e., p s ~ p. 

Both, the QCD sum rule and the chiral EFT fields are well comparable in terms of 
a density expansion since both are obtained to leading order in density. In the case of 
the sum rules this corresponds to a Fermi gas of non-interacting nucleons. To go beyond 
the Fermi gas approximation would require to include higher order terms in the operator 
product expansion and the density expansion of the condensates Q, |^ . In the EFT case 
higher orders in density can be introduced by a self-consistent dressing of the interaction 
(see discussion in Sec. IVI|) and of course by higher order terms in perturbation series which 
would finally end up in a full resummation of the Brueckner ladder diagrams. 

At moderate nuclear densities the agreement between the QCD sum rules and N 3 LO is 
quite remarkable. At higher densities the results from the sum rules tend to overshoot the 
N 3 LO values which is, however, not too astonishing since the relations (J51|) are valid in the 
low density limit. 
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FIG. 11: (Color online) Density dependence of the tree level scalar and vector self-energy compo- 
nents in nuclear matter obtained with the chiral EFT NN interaction 2^|. The fields obtained 
in leading order (LO) up to next-to-next-to-next-to-leading order (N 3 LO) are shown. The results 
from leading order QCD sum rules are shown as well. 



In view of the fact that in chiral NN dynamics the fields are dominantly generated by 
NLO contact terms, one could be tempted to interpret the present results in the way that 
the reduction of the quark condensates occurs at NLO in the chiral expansion. However, 
as discussed above such an interpretation is not straightforward. A closer inspection of the 
terms which drive the sum rule result reveals the following: the coefficient A in (|54|) is related 
to the unknown coupling C\ in the effective ChPT pion-nucleon Lagrangian (48 1. Becher and 
Leutwyler extracted a value of A = 3.7 GeV -1 fitting the elastic nN scattering amplitude at 
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threshold |49(. Inserting this value into the sum rule expression corresponds to a scalar 
self-energy (at kp = 1,35 fm _1 ) of S s = —513 MeV at order mj. At order m f , i.e., when 
the 0(m%) term in the expansion (J54)) is included, the sigma term of 46.7 MeV is already 
close to its empirical value and a self-energy of S s = —340 MeV is obtained. Although this 
value for E s is astonishingly close to the NLO result from chiral NN scattering, one has to 
keep in mind that already the LO result is of order m% in the pion mass. In contrast to the 
sum-rule approach there appears no significant repulsive contribution from pion dynamics 
which would correspond to the 0(m\) term in (|54j) . 



The present results are therefore in qualitative agreement with the findings of Refs. 
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421 ]. namely that long-distance physics related to pion dynamics plays only a minor role for 



the reduction of the nucleon mass in matter. Relating the in-medium nucleon mass to the in- 
medium scalar condensate through expression (J53*Jl one should be very careful. Although the 
sum rule mean fields, Eqs. (J49|) and ([50)1 . provide a reasonable approximation to the mean 
fields from chiral EFT, both approaches do not reflect the same physical concepts. The 
sum rule approach assumes that the nucleon properties are determined by the interaction 
with the in-medium condensates while conventional many-body approaches assume that the 
in-medium properties are determined by the interaction between the nucleons. 



VI. EQUATION OF STATE 

Until now all calculations in this paper have been performed at tree level. It is, however, 
a well known fact that a realistic description of nuclear dynamics requires correlations be- 
yond Hartree-Fock. Short-range correlations are known to be essential for nuclear binding 
whenever realistic interactions are used. This leads in lowest order of the Brueckner hole-line 
expansion to the ladder approximation of the Bethe-Goldstone equation for the in-medium 
G-matrix Q], or the Bethe-Salpeter equation in the relativistic case jljj. In contrast to 
non-relativistic BHF where the saturation points of isospin saturated matter are allocated 
on the so-called Coester line, the relativistic Dirac-Brueckner-Hartree-Fock approach leads 



to rather reasonable saturation properties 



12 



13. 



14 1 . For a review see |5 



In Hartree-Fock the matter turns out to be unbound, in particular when high precision 
potentials with a relatively strong repulsive hard core are applied, e.g. OBE type potentials 
or Argonne Vi$- The situation is qualitatively different for low momentum interactions 
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(Mow k, Idaho N 3 LO) where the hard core is strongly suppressed by the high momentum 
cut-offs. For these interactions isospin saturated nuclear matter collapses and Brueckner 
ladder correlations do not improve on this situation j^J. Here the matter has to be stabilised 
by the inclusion of repulsive three-body-forces j5l|. Doing so, there appears a strong cut-off 
dependence at tree-level which can be removed when the second order term of the Brueckner 
perturbation series is added. V\ ow k in combination with three-body-forces does not require 
a full resummation of the ladder diagrams but can already be treated within second-order 
perturbation theory jsi|. 

In the present work we do not aim for a fully realistic description of the nuclear many- 
body problem but restrict the investigations to the Hartree-Fock level. The tree-level results 
discussed up to now are of leading order in density p. Higher order corrections in density can 
be taken into account when the bare potential matrix elements are replaced by in-medium 
matrix elements V i— > V*. In the relativistic approach such a treatment is well defined. It 
means to evaluate the corresponding Feynman amplitudes (jSJ) through an in-medium spinor 
basis u* x {k) where the nucleons are dressed by the self-energy. Such a treatment requires, 
however, a definite structure of the interaction which allows to evaluate corresponding in- 
medium amplitudes. It is therefore at present restricted to OBE-type potentials. 

The dressing of the interaction through the self-energy leads automatically to a self- 
consistency problem which is e.g. solved within DBHF. The higher order density dependences 
which are introduced by such a procedure are considered to be one of the essential reasons 
for the improved saturation behaviour of relativistic DBHF compared to non-relativistic 
BHF. In the following we will study the role of self-consistency at the Hartree-Fock level. 

As already mentioned, in a relativistic framework one uses an in-medium spinor basis 
where the scalar and vector self-energy components from Sec. IIVI enter via effective masses 
and momenta, see Eq. ([46)1 . Furthermore the spatial vector self-energy component is usually 
absorbed introducing reduced effective masses and momenta 

M* ~ k* 
M* = , k* = ^— . (56) 

Hence the kinetic energy can be written as 



k* = E* = — ^— = \/k 2 + M* 2 (57) 
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and the in-medium spinors of helicity A are given by 



ul(k) 



E* + M 
\l 2M* 




(58) 



Thus the effective mass M* introduces a density dependence into the interaction. The effec- 
tive mass is, however, in general not only density but also momentum dependent. Based on 
the observation that this explicit momentum dependence is moderate, it is usually neglected 
and M* is fixed at the reference point |k| = kp. In the so-called reference spectrum approxi- 
mation the reduced effective mass M F = M*(|k| = kp, kp) serves as an iteration parameter. 
M* is then the solution of the non-linear equation 



which follows from the formulae above. Self-consistency is now achieved by determining 
for a given start value of M* the in-medium matrix elements V^^q^q). Therefore the 
Lorentz invariant amplitudes F^(|q|,6>) and g^dq], 9), Eqs. (J29"j) and (J3~2j) . as well as the 
transformation matrix Cj m of Eq. ()31|) depend on M* and the Fermi momentum kp since 
the plane- wave helicity states |AiA2q) of Eq. (fHU|) are now medium- dependent (f5%|) . The 
next step is to compute the self-energy components S s , S and kS v . Since the Dirac prop- 
agator (}37))) describes dressed quasi-particles now, also in (jUJ, (}4*2|) and (f4*Hj) the mass M 
and energy E have to be replaced by the effective quantities M* , E*. Finally the new M* 
is determined. This iteration procedure is repeated until convergence is reached. 

In Fig. El the results for the self-consistently calculated self-energy components S s and 
Eo for Bonn A, Nijm 93 and Nijmegen I are shown as a function of the Fermi momentum 
and compared to the tree level results from Fig. El Fo r the Bonn A case the result of a full 
self-consistent DBHF calculation is shown as well jl4] . From this figure two features can be 
observed: the higher order density dependences which are introduced by the dressing of the 
potential affect mainly the scalar part of the self-energy. The modifications of S are moder- 
ate while S s is significantly reduced. The short range ladder correlations included in the full 
DBHF calculation (Bonn A) influence the self-energy in an opposite way. The deviations of 
E s from the self-consistent HF result are rather small, however, the vector component gets 
now strongly suppressed. This fact is understandable since the ladder correlations prevent 
the two-nucleon wave functions from too strong overlap with the hard core. In OBE poten- 



M* = M + S s (fc F , M*) - M*S v (fc F , M*) 



(59) 
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FIG. 12: (Color online) Comparison of the tree level scalar and vector self-energy components 
(dashed line) with self-consistent results (solid line). Additionally a full self-consistent DBHF 
calculation is shown in the first graph denoted by dots. 



tials the hard core is mainly mediated by vector w-exchange and determines thus the vector 
self-energy component. 

With the self-consistent Hartree-Fock self-energies at hand one can now determine the 
equation of state (EOS). Like in DBHF the EOS, i.e., the energy per particle is defined as 
the kinetic plus half of the potential energy 



1 1 M* 

E/A=- £<^( k )l7 .k + M + -£(A;)K(k)>-^-M 

P k.A 2 E 



p Jf 2tt 3 



((1 + E v (|k|))£*-E°(|k|)) 



2E* 



[E s {\k\)M* - E M (|k|)£*" 



(60) 
M (61) 



with the self-consistent spinors u\ from Eq. (|58p. 

In Fig. El we present the self-consistent Hartree-Fock results for the energy per particle 
in symmetric nuclear matter calculated from the Bonn A, Nijm93, Nijmegen I potentials as 
a function of the Fermi momentum kp which is a measure for the density p = 2/(3n 2 )k F . 

Also a non-self-consistent calculation is shown (dashed line) where the energy per particle 
is given by 



E/A 



1 r d 3 k 

p Jf 2-7T 3 



(62) 



with U B _ p _(k, kp) as defined in Eq. (|43|). In this case one obtains the same result as in a non- 
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FIG. 13: (Color online) Hartree-Fock calculation of the nuclear equation of state, i.e., energy per 
particle E/A as a function of the Fermi momentum hp for three different potentials. The dashed 
line indicates a tree level calculation and the solid line represents a self-consistent Hartree-Fock 
calculation, i.e., higher order corrections in density are included. 



relativistic Hartree-Fock calculation (denoted by stars in Fig. fl~3|) . The latter demonstrates 
again the numerical accuracy of the procedures. 

For the Bonn A case again the EOS from the full DBHF calculations is shown as a 
reference Q. It is dear that .adder correlations and other in-medinnr effects such as Pan.,- 
blocking of intermediate states in the Bethe-Salpeter equation are responsible for nuclear 
saturation. The relatively moderate deviations from self-consistent Hartree-Fock at the 
scale of the self-energies in Fig. El are essential at the scale of the binding energy. Like 
in relativistic mean field theory of QHD subtle cancellation effects in the large scalar and 
vector fields are responsible for nuclear binding. 

The higher order density dependences introduced via the dressing of the bare interaction 
V lead to significantly more repulsion at the level of the EOS. This is a direct consequence of 
the reduced attractive scalar field (see Fig. I12J). Thus Fig. Unserves also as a demonstration 
for the success of DBHF compared to BHF what concerns the quantitative description of nu- 
clear saturation: In particular for modern high precision potentials such as Bonn, Nijmegen 
or Argonne v\g the BHF approach leads to strong over-binding and too high saturation 
densities. The additional repulsion introduced by higher order terms in density through 
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the dressed potentials shifts the corresponding saturation points towards the empirical re- 



gion 



13, 



16, 55]. to stress that the density dependence of the dressed potential V* 

should not be mixed up with the density dependence of the G-matrix. The latter originates 
from the dressed two-nucleon propagator and the Pauli-operator in the Bethe-Goldstone (or 
Bethe-Salpeter) equations while V* enters into the Bethe-Salpeter for iteration. In non- 
relativistic BHF or variational calculations j3, El a non-linear density dependence which 
improves the saturation behaviour is usually introduced through net repulsive three-body- 
forces. In such a treatment the dependence on the third particle is integrated out such that 
one is left with an additional effective density dependent two-body force which acts in a 
similar way as a dressing of the two-body interaction. In this context one should mention 
that a dressing of the interaction has also more subtle consequences when iterated in the 
Bethe-Salpeter equation. It leads e.g. to a quenching of the second order OPE exchange 5^ 
which plays an essential role for saturation in non- relativist ic approaches. 

In summary, one could expect that a dressing of the interaction would allow to comply 
with weaker three-body forces which may in particular be of interest concerning the appli- 
cation of low momentum EFT potentials to the nuclear many-body problem. As the studies 



of Bogner et al. [51| have demonstrated, V\ ow k requires rather strong three-body forces in 
order to stabilise nuclear matter. There the strength of the three-body contributions has 
already been pushed to its limits. Although a dressing of the interaction will probably not 
be possible for Vi ow k due to the partially non-analytic structure of the potential, it may be 
a promising perspective for the application of other EFT potentials, e.g. the chiral N 3 LO. 



VII. SUMMARY 



The appearance of large scalar and vector fields is a well established feature of relativistic 
nuclear dynamics. The saturation mechanism of nuclear matter or the single particle poten- 
tial in finite nuclei are obtained by subtle cancellation effects between large attractive scalar 
and repulsive vector fields. These fields occur already at tree level and do not change too 
much when realistic many-body calculations are performed. Full self-consistent Brueckner 
calculations which account for short-range ladder correlations lead to mean fields of similar 
size, i.e., of several hundred MeV magnitude. The size of the scalar- vector fields coincides 
with the values derived from relativistic mean field phenomenology by fits to finite nuclei. 



43 



Alternatively, QCD sum rules come to the same results. 

The present work addresses the question about the origin of these fields. When the 
nucleon-nucleon interaction is described within the framework of a meson exchange picture, 
the situation is rather clear. The Lorentz character of the mesons determines automatically 
the Lorentz character of the interaction at the corresponding scale: the short-range repulsion 
is due to vector exchange (u, p) while the intermediate range attraction originates from scalar 
exchange (a). As a direct consequence this leads to the existence of large scalar and vector 
mean fields in nuclear matter. However, these fields are not observables. It is therefore 
a fundamental questions of nuclear physics whether the appearance of large scalar/vector 
fields is intimately connected to the meson exchange picture or if it is a general consequence 
of the vacuum NN interaction. 

To address the question in a model independent way, we based the present study on a 
broad set of modern high precision NN potentials: Bonn, CD-Bonn, Nijmegen, Argonne 
t>i 8 , Reid93, Idaho N 3 LO and V\ ow k- Except the fact that all these potentials fit NN 
scattering data with high accuracy they are based partially on quite different theoretical 
concepts, i.e., the traditional meson exchange picture (Bonn, CD-Bonn, Nijmegen), a purely 
phenomenological philosophy (Argonne fig, Reid93) or QCD inspired effective field theory 
approaches (Idaho N 3 LO, Vj ow k)- 

For this purpose the potentials were projected on a relativistic operator basis in Dirac 
space. This was achieved using standard projection techniques which transform from a 
partial wave basis, i.e., the basis where the potentials are originally given, to the basis of 
covariant amplitudes in Dirac space. The idea behind this approach is that both, relativistic 
and non-relativistic descriptions of the NN interaction have common features, i.e., they are 
based on a certain operator structure in spin-isospin space and invoke certain scales: the 
long-range part of scale m n , essentially given by one-pion exchange, the intermediate range 
attraction and the short-range repulsion. In the meson exchange picture the various scales 
are associated with the meson masses which mediate the interaction. The various approaches 
can now be compared at the level of these covariant amplitudes where we observe a remark- 
able agreement between the meson exchange potentials (Bonn, CD-Bonn, Nijmegen), the 
phenomenological non-relativistic potentials (Argonne i>i 8 , Reid93) and the EFT potentials 
(Idaho N 3 LO, V lowk ). 

Moreover, this procedure allows now to calculate the relativistic self-energy operator in 
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nuclear matter. The key result of the present investigations is the tree level self-energy in 
nuclear matter. The structure of the nucleon-nucleon interaction enforces the existence of 
large scalar and vector fields. This is found to be a model independent fact, true for all 
types of interactions which have been considered. The scale of these fields is set at tree 
level. Although essential for nuclear binding and saturation, higher order correlations, in 
particular short-range correlations, change the size of the fields by less than 25%. The 
magnitude of the tree-level fields is very similar to that predicted by relativistic mean field 
phenomenology and relativistic many body calculations. 

The connection to QCD as the underlying theory of strong interactions is established 
by chiral effective theory. EFT nucleon-nucleon potentials are derived from a systematic 
expansion of an effective Lagrangian which respects the basic symmetries of QCD. Chiral 
EFT is considered as the exact mapping of QCD on effective hadronic degrees of freedom 
in the non-perturbative regime. Subjecting the chiral N 3 LO Idaho potential to the present 
projection scheme we can make the following statements: In nuclear matter scalar and 
vector mean fields of the same sign and magnitude are generated as by the meson exchange 
or phenomenological potentials. These fields are generated by contact terms which occur 
at next-to-leading order in the chiral expansion. These are four-nucleon contact terms with 
two derivatives which generate the short-range spin-orbit interaction. The strength of the 
corresponding low energy constants, in particular those connected to the spin-orbit force, is 
dictated by P-wave NN scattering data. Pion dynamics as well as LO and N 3 LO contacts 
provide only corrections to the fields generated by the NLO contact terms. EFT is therefore 
in perfect agreement with Dirac phenomenology where it is known since a long time that the 
large scalar/vector fields are generated by the short-range vector (u) and scalar (a) mesons 
which are connected intimately to the large spin-orbit interaction. We conclude that this is 
a direct consequence of P-wave NN scattering. 

Like in OBE models and RMF theory, in EFT the reduction of the nucleon mass M* = 
M + S s is driven by short- distance physics. Long-distance physics from virtual pions, i.e. 
the non- analytic term in the expansion of <jn gives a sizable contribution to the modification 
of the in-medium quark condensate. Such contributions are, however, found to play only a 
minor role for the reduction of the nucleon mass. Nevertheless, at moderate nuclear densities 
the N 3 LO scalar and vector fields agree almost perfectly with the prediction from leading 
order QCD sum rules. For future perspectives chiral EFT in combination with projection 
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techniques may allow to determine the relativistic anti-proton potential in matter in a model 
independent way. Here the meson-exchange picture predicts a change in sign of the vector 
field due to g-parity and hence an extremely deep attractive potential. Such investigations 
in particular will be interesting in view of the forthcoming anti-proton facilities, e.g. Panda 
at FAIR |59|. 

Finally we investigated implications of higher order corrections in density on the nuclear 
EOS. A dressing of the potential through self-consistently determined self-energies leads 
already at the Hartree-Fock level to significantly more repulsion in the EOS. At present 
these investigations were restricted to OBE type potentials. But to include such higher 
order terms in density might open a promising perspective also for EFT potentials when 
applied to the nuclear many-body problem. 
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